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Abstract 

Neutral pion (tt^) production is calculated in a leading order (LO) perturbative 
QCD-based model in pp and dAu collisions at ^/s = 200 AGeV at midrapidity. 
The model includes the transverse component of the initial parton distribution. We 
compare our results for pp collision to experimental data at RHIC energy. We repeat 
our calculation for the dAu collision and investigate the interplay between shadowing 
and multiple scattering. In central dAu collisions the influence of possible jet energy 
loss on cold nuclear matter is discussed and numerical results are displayed. 

Key words: Nuclear multiscattering, nuclear shadowing, energy loss 
PACS: 12.38.Bx, 13.87.-a, 24.85. +p, 25.75.-q, 25.75. Gz 



1 Introduction 



Experimental data on 7r° production in Au + Au collisions at ^/s = 200 AGeV 
display a very strong suppression pattern in the transverse momentum region 
2 GeV < pt < 20 GeV [1], which indicates a clear evidence of the presence of 
induced jet energy loss in hot dense matter [2-5]. Quantitative calculations 
of the energy loss in Au + Au collisions require the precise description of 
the available data from d + Au collisions. These data incorporate initial state 
effects (nuclear multiscattering and nuclear shadowing) and possible final state 
effects (e.g. energy loss in cold nuclear matter), which should be understood. 
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2 The parton model for pp collisions 

In this paper we accomplish the analysis of pp and dAu data in a leading-order 
(LO) perturbative QCD based model [6,7]. We plan to repeat our analysis in 
next-to-leading (NLO) order [8], when final data will be available for pion 
production in dAu collisions. The invariant cross section for pion production 
in pp collision can be described in a LO pQCD-improved parton model on the 
basis of the factorization theorem as a convolution [6,7]: 




where da/di is the hard scattering cross section of the partonic subprocess 
ab cd, and the fragmentation function (FF), Dh/c{zc,Q''^) gives the proba- 
bility for parton c to fragment into hadron h with momentum fraction Zc at 
scale Q'. For fragmentation functions we use the KKP parameterization [9]. 

The functions fa/p{x, k^a, Q^) and fb/pix, k^b, Q^) are the parton distribution 
functions (PDFs) for the colliding partons a and b in the interacting protons 
as functions of momentum fraction x, at scale Q. We use a product approx- 
imation, namely /(x, kj-, Q^) = f{x,Q'^)g{'kT), where the function f{x,Q^) 
represents the standard LO PDF (here it is the GRV set [10]). In our phe- 
nomenological approach the transverse-momentum distribution is described 
by a Gaussian, 

Here, (k^) is the 2-dimensional width of the kx distribution and it is related 
to the magnitude of the average transverse momentum of a parton as {k^) — 

The validity of this approximation and the introduction of the transverse 
component into the PDFs were under debate for a long time. However, recent 
measurement of di- hadron correlations by PHENIX [11] and the theoretical 
investigation of these data [12] have shown that a large value for the kx- 
imbalance can be extracted from existing di-hadron data. The origin of this 
/cr-imbalance is three-fold: intrinsic transverse momenta of the partons, non- 
perturbativc soft radiation and higher order radiation terms. We consider these 
effects phenomenologically and our (A;|,) value accumulates all of these effects. 
As we will see, the size of this (k^) can be extracted consistently from one- 
particle spectra and two-particle correlations. 
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In our present study we consider fixed scales: the factorization scale is con- 
nected to the momentum of the intermediate jet, Q = hz-pg (where Pq = pt/z^.), 
while the fragmentation scale is connected to the final hadron momentum, 
Qf = K ■ pt- The value of n can be varied in a wide range, however we will 
use K = 2/3, which is the best choice to obtain agreement between data and 
theory at high-py, where non-perturbative effects have a small contribution. 

Figure 1 displays our 
results at = 200 GeV 
for pp TT^X [13] at 

different values of the 
fc-r— imbalance. Without 
transverse momentum 
component for partons 
the theoretical results 
underestimate the data by 
a factor of 2 in the window 
2 GeV < Pt < 5 GeV. 
To minimize the difference 
between RHIC data and 
the model, a value of 
(k^) = 2.5 GeV2 can be se- 
lected. This value indicates 

iphpa^r) = m) = 5 GeV^ 
which is compatible with 
the value obtained in 
Fig. 1. (Color online) One-particle spectra for vr'^ pro- Refs. [11,12] from di- 
duction in pp collisions at ^/s = 200 GeV at different hadron spectra, 
intrinsic /cj" values compared to PHENIX data [13]. 

Our earlier studies have shown that the theoretical reproduction of the one- 
hadron spectra in a pQCD frame strongly requires the introduction of a large 
/ct— imbalance in both LO and NLO calculations [14]. The analysis of di- 
hadron correlations in NLO frame would give more insight into this problem. 



3 The parton model for dAu collisions 

Considering the dAu collision, the hard pion production cross section can be 
written as an integral over impact parameter b, where the geometry of the 
collision is described in the Glauber picture [6,7]: 

E,^-fl . /d'.dV Ur) .,„(|b - .1) . E, ^'^^-t- ■ (3) 

a p J dp 
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where the proton-proton cross section on the right hand side represents the 
cross section from eq. (1), but with the broadened widths of the transverse- 
momentum distributions in eq. (2), as a consequence of nuclear multiscattering 
(see eq. (4)). In eq. (3), tAu{b) = J dz pAu{b, z) is the nuclear thickness function 
(in terms of the density distribution of the gold nucleus, pau), normalized as 
/ (^^htAuip) = ^Au = 197. For the deuteron, one could use a superposition 
of a pAu and a nAu collision, or a distribution for the nucleons inside the 
deuteron. Here we apply a hard-sphere approximation for the deuteron with 
A— 2 for estimating the nuclear effects. Also, since 7r° production is not sensi- 
tive to isospin, we continue to use the notation "pA" when talking about the 
interaction of any nucleon with a nucleus. 

The initial state broadening of the incoming parton distribution function is 
accounted for by an increase in the width of the gaussian parton transverse 
momentum distribution in eq. (2) [7,15]: 

(4>pA = {kl)pp + C ■ hpAib) . (4) 



Here, {k^)pp is the width of the transverse momentum distribution of partons 
in pp collisions, hpA(h) describes the number of effective nucleon- nucleon (NN) 
collisions at impact parameter 6, which impart an average transverse momen- 
tum squared C. The effectivity function hpA{b) can be written in terms of the 
number of collisions suffered by the incoming proton in the target nucleus, 
hpA{b) — i'A{b) — 1. Here, UA^b) = aMNtAip), with a^M being the inelastic 
nucleon-nucleon cross section. For the factor C and z^a(^) we will use the find- 
ings of Ref . [7] , where the systematic analysis of pA reactions was performed 
in LO and the characteristics of the Cronin effect were determined at LO 
level. Following Ref. [7], we assume that only a limited number of semi-hard 
collisions (with maximum i'A{b)max — 4) contributes to the broadening, and 
we have found C = 0.4 GeV^. The maximum of the broadening is ~1 GeV^ 
and this value determines the maximum increase in the particle yield. 

It is well-known that the PDFs are modified in the nuclear environment and 
nuclear shadowing appears for partons with momentum fraction x < 0.1. This 
is taken into account by various shadowing parameterizations [16-19]. In the 
present work, we display results obtained with the EKS parameterization, 
which has an antishadowing feature[16], and with the updated HIJING pa- 
rameterization[17], which incorporates different quark and gluon shadowing 
and requires the introduction of nuclear multiscattering. As a third case, we 
consider the nuclear parton distribution functions introduced by HKM [18]. 
Since beyond nuclear shadowing the nuclear multiscattering may appear, we 
will apply the HKM parametrization alone and in a combination with our 
multiscattering description. 
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4 Results on pion production in dAu collisions 



Including the multiscattering and shadowing effects one can calculate the in- 
variant cross section for pion production in d + Au collision. Moreover, intro- 
ducing the nuclear modification factor RdAu, as 



_ £'^dcr^^"/d^p _ £'^d(T^^''( with nuclear effects)/d^p 
'^^'^ ~ Nun ■ EAal^/d^P ~ ^^da^^«( no nuclear ef f ects)/d^p ^ 

nuclear effects can be investigated clearly and efficiently, using a linear scale. 
The value of can be determined using the Glauber geometrical overlap 
integral as in eq. (3). Here we apply the right-most equation in our theoret- 
ical calculations, which does not require the determination of N^n from the 
Glauber model. 
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Fig. 2. (Color online) Nuclear modification factor for -k^ production in dAu collisions 
at ^/snn = 200 AGeV at different centralities. Theoretical results (see text for 
details) are compared to preliminary PHENIX data [20]. 

Figure 2 summarizes our theoretical results in four centrality bins for the nu- 
clear modification factor RdAu in dAu collisions at ^/s = 200 AGeV, displaying 
different shadowing models. Preliminary PHENIX data are from Ref. [20]. 
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At first glance, the experimental data have a striking message at all central- 
ity bins: the nuclear modification factor shows a ±10% deviation from unity 
in wide momentum regions in all centrality bins, including the most central 
and the most peripheral ones. Since these ±10% deviations overlap with the 
systematic error of these experimental data, we can not exclude the scenario 
of a null effect, which means no multiscattering and no shadowing are present 
at RHIC energies. Fortunately the nuclear modification factor has been mea- 
sured at higher rapidities and the influence of nuclear shadowing can be seen 
clearly [21], which excludes the null effect scenario. This means that there is 
a delicate interplay between nuclear shadowing and multiscattering (or the 
antishadowing feature of the shadowing functions) . 

Investigating the experimental data from this point of view, we may recognize 
a slight b-depcndcnce for both multiscattering and shadowing. In the most 
peripheral collisions (60-88 %) multiscattering overwhelms the shadowing - 
or shadowing function has a clear antishadowing region to be involved in 
the pion production in the region 3 < pr < 7 GeV/c. The enhancement 
can be reproduced by the EKS (thin dashed line) and the HKM shadowing 
(solid line), the scenario of 'HKM+multiscattering' (dotted line) is satisfactory 
only because of the large error bars. HIJING (thick dashed line) displays no 
nuclear modification — more precise data may help us to choose between the 
different scenarios. In the most central collisions (0-20 %) in the intermediate 
transverse momentum window the +10% enhancement appears through the 
superposition of a stronger shadowing and a more effective multiscattering 
connected to the geometry of the collisions. If final data will have a smaller 
error bars, then we may be able to select between the different theoretical 
models. 

However, data in the most central bins (0-20 %) display an interesting phe- 
nomenon at high-p^: although all model calculations with nuclear multiscat- 
tering and shadowing lead to unity, the experimental data are decreasing with 
increasing transverse momentum. The authors are not aware of any shadow- 
ing parametrization, which could yield such a suppression around = 6 — 10 
GeV at RHIC energies, and any nuclear multiscattering should work in the 
opposite direction. This suppression in central dAu collisions may indicate 
the presence of a new phenomena, namely a moderate jet energy loss in 
the cold matter produced in dAu collisions. Cold quenching of high energy 
quarks is known from the HERMES experiment [22], and recently it was in- 
vestigated theoretically (sec Refs. [23-25]). The question of cold quenching is 
very complicated, and it is strongly related to the microscopical description of 
hadronization (see Ref. [25]). Thus recent dAu data with large error bars can 
not give a deep insight into this phenomena beyond indicating its existence 
in the most central bins. We estimated this effect with a simple rescaling of 
the 'HKM+multiscattering' scenario. Figure 3 indicates that dAu data in the 
most central bin can accommodate a 10-15 % suppression. 
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Fig. 3. (Color online) Nuclear modification factor for vr'' production in central dAu 
collisions at y/sJm = 200 AGeV. Theoretical result (dotted line) is obtained by 
simple rescaling of the 'HKM+multiscattering' scenario (see text for details). Pre- 
liminary PHENIX data are from Ref. [20]. 



5 Summary 



We have investigated the interplay between nuclear shadowing and multi- 
scattering at RHIC energies in dAu collisions. Since the nuclear modification 
factor, RdAui has a value close to unity in large transverse momentum regions, 
the two effects seem to balance each other within ±15%, which makes the 
separation and the quantitative analysis of these effects in the mid-rapidity 
difficult. Final data with smaller error bars may lessen this complication. The 
characteristic suppression of the nuclear modification factor at high-pT- in the 
most central dAu collisions may indicate the presence of a moderate jet en- 
ergy loss. Quantitative analysis requires the introduction of proper theoretical 
models for jet energy loss in cold nuclear matter. 
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